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This paper generalizes and expands upon the work 1] where we introduced a scheme for fault- 
tolerant holonomic quantum computation (HQC) on stabilizer codes. HQC is an all-geometric 
strategy based on non-Abelian adiabatic holonomies, which is known to be robust against various 
, types of errors in the control parameters. The scheme we present shows that HQC is a scalable 

■ method of computation and opens the possibility for combining the benefits of error correction with 

' the inherent resilience of the holonomic approach. We show that with the Bacon-Shor code the 

scheme can be implemented using Hamiltonian operators of weight 2 and 3. 

^3 : 

! I. INTRODUCTION 

o : 

CN . There are two main sources of errors in quantum computers — environment-induced decoherence and imperfect 
' control. It has been shown that if the errors of each type are sufficiently uncorrelated and their rates are below a 
, certain threshold, it is possible to implement reliably an arbitrarily long computational task with a modest resource 
rS i' overhead [2,, ^, j5„ ^, 7, 8, 9]. This result, known as the quantum accuracy threshold theorem, is based on the idea 
^ \ of quantum error correction (QEC) [l^ [ll| — a universal software strategy to combat noise in quantum computers. 
' In addition to the software approach, there have also been proposals to deal with the effects of noise by hardware 
\ methods that provide robustness through their inherent properties. One such method is holonomic quantum compu- 
tation f HQC) an adiabatic, all-geometric method of computation that uses non-Abelian generalizations |1J| 

of the Berry phase [15|. It has been shown that due to its geometric nature, this approach is robust against various 

types of errors in the control parameters driving the evolution [l^, [13, [H, [H, [lOI, and thus provides a degree of 
• built-in resilience at the hardware level. 
^ , In Rcf. j21|| HQC was combined with the method of decoherence- free subspaces (DFSs) [13, [1^ [13], which was 
£]0 ■ the first step towards systematic error protection in conjunction with the holonomic approach. DFSs provide passive 
I protection, in the sense that no syndrome measurement and feedback are required. The DPS approach is applicable 
^vq I provided there is a symmetry in the system-bath coupling, which results in a sufhciently large decoupled subspace 
in which quantum information can be stored. Collective decoherence is a well-known example of such a symmetry 
(for a review see Ref [1^). In the absence of symmetries, the addition of syndrome measurements and feedback is 
useful and leads to the approach based on active error correction [l^, [ll| . Active error correction is also the basis of 
, quantum fault tolerance, which is essential for the scalability of any method of computation. Even if we assume that 
IL* ■ the system is perfectly protected from environment-induced errors, when the size of the circuit increases, errors due 
to imperfect operations would accumulate detrimentally unless they were corrected. Therefore, scalability of HQC 
requires combining the holonomic approach with active error correction. 
?H In Ref. (Ij, we reported a scheme which combines HQC with the techniques for fault-tolerant computation on 

stabilizer codes. This demonstrated for the first time that HQC is a scalable method of computation. The scheme 
uses Hamiltonians which are elements of the stabilizer — or in the case of subsystem codes, elements of the stabilizer and 
gauge groups. Encoded gates are implemented by slowly varying the Hamiltonians along suitable paths in parameter 
space, such that from the point of view of the basis of the full Hilbert space, the states in each eigenspace undergo the 
same transversal transformation. On certain codes such as the 9-qubit Shor code [lo| or its subsystem generalizations 
[26l. [27|. universal computation according to our scheme can be implemented with Hamiltonians of weight 2 and 3. 

This paper generalizes and expands upon the work We provide details on various proofs sketched in Ref. [l|, 
clarify many points and analyze properties of the scheme that were not discussed there. We examine in detail the 
construction for the Bacon-Shor code [H, [13|i discuss the adiabatic approximation for different parametrizations of 
the Hamiltonians, and provide explicit calculations of the holonomy in the implementation of the Z (phase-flip) gate 
for two different types of interpolation. 
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II. PRELIMINARIES 



A. Holonomic quantum computation 

Let {H\} be a family of Hamiltonians on an A^-dimensional Hilbert space, which is continuously parametrized by 
a point A in a control-parameter manifold M. Assume that the family has the same degeneracy structure, i.e., there 
are no level crossings. The Hamiltonians can then be written as H\ = X]n=i ^nW^nW: where {£nW}n=i the 
R different d„-fold degenerate eigenvalues of Hx, (X]n=i ~ -^)' ^^"-^ n„(A) are the projectors on the corresponding 
eigenspaces. If the parameter A is changed adiabatically, a state which initially belongs to an eigenspace of the 
Hamiltonian will remain in the corresponding eigenspace as the Hamiltonian evolves. The unitary evolution that 
results from the action of the Hamiltonian H{t) :— Hx(t) is 

C/(t) = rexp(-i rdri/(r)) =©f^ie-"Wt/l(i), (1) 



where w„(t) = — (iTe„(A(T)) is a dynamical phase, and the intrinsically geometric operators U^^{t) are given by 
the following path-ordered exponents: 

rHt) 

U^Jt)=Vexp{ An). (2) 
Jx{o) 

Here An is the Wilczek-Zee connection 'T^ for the n-th eigenspace. An = ^fj, An,fj.dX^ , where ^n,/^ has matrix 
elements [l^ 



d 

{An,f,)ap = {na;X\-^\nf3;X). (3) 



The parameters A^ are local coordinates on dimAl) and {|na;A)}^^]^ is an orthonormal basis of the 

n-th eigenspace of the Hamiltonian at the point A. 

When the path A(t) forms a loop 7(t), 7(0) = jiT) ~ Aq, the unitary matrix with components {U2)ap appearing in 



^iU:i)o,p\na;O){n0;O\=U^jT)^Vexpi(f An) (4) 

Q/3 "'t 



is called the holonomy associated with the loop. In the case when the n-th energy level is non-degenerate (c?„ = 1), the 
corresponding holonomy reduces to the Berry phase [Tsj . The holonomy is a geometric object which is invariant 
under gauge transformations corresponding to changing the basis of the n-th eigenspace along the curve 7. The set 
Hol(A„) = {t7T;|7 e LxoiM)}, where Lx^iM) = {7 : [0, T] ^ M\-fiO) = -/(T) = Aq} is the space of ah loops based 
on Ao, is a subgroup of U{dn) called the holonomy group. 

In Refs. [13, [3] it was shown that if the dimension of the control manifold is sufficiently large, quantum holonomies 
can be used as a means of universal quantum computation. In this approach, logical states are encoded in the 
degenerate eigenspace of a Hamiltonian and gates are implemented by adiabatically varying the Hamiltonian along 
suitable loops in the parameter manifold (for a construction of a universal set of gates, see also Ref. [28]). 

We point out that many assumptions behind this simple model of HQC can be relaxed. For example, if we are 
interested in performing computation in the n-th eigenspace of the Hamiltonian, it is sufficient that this eigenspace 
is adiabatically decoupled from the rest of the Hilbert space, and it is not necessary that there are no crossings 
between other energy levels that are separated from the n-th level by energy gaps. Furthermore, in order to obtain a 
gauge- invariant expression for the geometric transformation taking place inside the n-th eigenspace, it is not necessary 
that the entire Hamiltonian undergoes a cyclic change — it is enough to take the n-th eigenspace around a loop. In 
fact, the form of the restriction of the Hamiltonian on the orthogonal complement of that subspace is irrelevant since 
what is important for the geometric transformation taking place inside an adiabatically decoupled eigenspace is how 
this subspace changes inside the full Hilbert space. More precisely, adiabatic quantum holonomies inside the n-th 
eigenspace can be equivalently understood as arising from parallel transport of vectors along the tautological bundle 
whose base is the Grassmannian parametrizing the set of d„-dimensional subspaces of the full Hilbert space, rather 
than from parallel transport along the corresponding bundle over the full space of control parameters. We note that 
even the requirement for a closed loop in the Grassmannian can be relaxed using the notion of open-path holonomies 
[29! . The approach that we pursue in this paper can be best understood as based on closed loops in the Grassmannian, 
even though — with a small modification — it can be made to be exactly of the type discussed in the above formulation, 
where the Hamiltonian family has a fixed degeneracy structure, and gates are implemented by loops in the control 
manifold (see Sec lII Cp . 
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B. Stabilizer codes and fault-tolerant computation 



33j. 



A large class of quantum error-correcting codes can be described by the so called stabilizer formalism [30|, |31| 
A stabilizer S is an Abelian subgroup of the Pauli group Gn on n qubits that does not contain the element — / 
The Pauli group consists of all possible n-fold tensor products of the Pauli matrices cr^ = X , = Y , = Z 
together with the multiplicative factors ±1, ±i. The stabilizer code corresponding to S is the subspace of all states 
\ip) which are left invariant under the action of every operator in S (G'|'0) = VG S S). It is easy to see that the 
stabilizer of a code encoding k qubits into n has n — k generators. Recently, a more general notion of codes has been 
introduced — subsystem or operator codes [13, — that employs the most general encoding of information, encoding 
in subsystems [s^ [13, [ll]. In the case of a single subsystem, the Hilbert space decomposes as H = ® (B K., 
where Ti.^ is the subsystem in which the logical information is encoded, Ti.^ is the gauge subsystem, and K, is the 
rest of the Hilbert space. For operator stabilizer codes, the stabilizer leaves the subspace TL^ ® Ti.^ invariant but 
the encoded information is invariant also under operations that act on the gauge subsystem. An operator stabilizer 
code encoding k qubits into n with r gauge qubits has n — r — k stabilizer generators, while the gauge group has 2r 
generators [s^]. According to the error-correction condition for stabihzer codes [ll,[3i|, a set of errors {Ei} in Gn 
(which without loss of generality are assumed to be Hermitian) is correctable by the code if and only if, for all i and 
j, EiEj anticommutes with at least one element of S, or otherwise belongs to S or to the gauge group. In this paper 
we will be concerned with stabilizer codes for the correction of single-qubit errors and the techniques for fault-tolerant 
computation [3, [3, i, [1, [1, [3, [S H on such codes. 

A quantum information processing scheme is called fault-tolerant if a single error occurring during the implemen- 
tation of any given operation introduces at most one error per block of the code . This property has to apply for 
unitary gates as well as measurements, including those that constitute the error-correcting operations themselves. 
Fault-tolerant schemes for computation on stabilizer codes generally depend on the code being used — some codes, 
such as the Bacon-Shor subsystem codes [H, [l^l, for example, are better suited for fault-tolerant computation than 
others [13] ■ In spite of these differences, however, it has been shown that fault-tolerant information processing is 
ossible on any stabilizer code [3, [3- The general procedure can be described briefly as follows. DiVincenzo and Shor 
demonstrated how to perform fault-tolerant measurements of the stabilizer for any stabilizer code. Their method 
makes use of an approach introduced by Shor [3|, which involves the "cat" state (|0...0) -|- \l...l))/^/2 which can be 
prepared and verified with a satisfactory precision. As pointed out by Gottesman [^, by the same method one can 
measure any operator in the Pauli group. Since the encoded X, Y and Z operators belong to the Pauli group for any 
stabilizer code [^, one can prepare fault tolerantly various superpositions of the logical basis states |0) and |1), such 
as 1+) = (jo) -I- |l))/-\/2, for example. The latter can be used to implement fault tolerantly the encoded Phase and 
Hadamard gates as long as a fault-tolerant C-NOT gate is available |8| . Gottesman showed how the C-NOT gate can 
be implemented fault tolerantly by first applying a transversal operation on four encoded qubits and then measuring 
the encoded X operator on two of them. Finally, for universal computation one needs a gate outside of the Clifford 
group, e.g., the Toffoli gate. The Toffoli construction was demonstrated first by Shor in Ref. [2,] for a specific type of 
codes — those obtained from doubly-even self-dual classical codes by the Calderbank-Shor-Steane (CSS) construction 
[4l|,[42|. Gottesman showed [3| that a transversal implementation of the same procedure exists for any stabilizer code. 



C. Overview of the scheme 



Note that the described method for universal fault-tolerant computation on stabilizer codes uses almost exclusively 
transversal operations — these are operations for which each qubit in a block interacts only with the corresponding 
qubit from another block or from a special ancillary state such as Shor's "cat" state (see also Steane's [i^ and Knill's 
[4J] methods). However, transversal operations are not the most general class of operations that do not lead to 
propagation of errors. For example, every transversal operation in a given fault-tolerant protocol can be substituted 
by the same operation followed by an operation that multiplies each syndrome subspace by a different phase, and the 
resultant protocol will still be fault-tolerant. This can be easily seen from the fact that if after a transversal operation 
the state is correctable, then it will still be correctable after multiplying each syndrome subspace by a phase because 
the correction procedure involves a projection on one of the syndrome subspaces and thus the overall phase in that 
subspace is irrelevant. An operation which is equal to a transversal operation followed by a transformation on the 
gauge subsystem can be similarly seen to be fault-tolerant. It is these more general fault-tolerant transformations by 
means of which we will realize fault-tolerant HQC. 

To explain the main idea behind our approach, let us consider the case of standard (subspace) stabilizer codes first. 
Our goal will be to find a holonomic realization of a universal set of encoded gates by adiabatically transporting the 
code space along suitable loops via sequences of elementary fault-tolerant transformations of the above type. The 
scheme we will present can be roughly described as follows. We choose as a starting Hamiltonian an element of the 
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stabilizer and vary this Hamiltonian in an adiabatic manner along appropriate paths so that from the point of view 
of the basis of the full Hilbert space, the vectors in both eigenspaces of the Hamiltonian undergo the same transversal 
transformation. Under this procedure, each eigenspace will acquire a dynamical phase corresponding to the energy 
of that eigenspace, but since the Hamiltonian is an element of the stabilizer, these phases will amount to relative 
phases between different syndrome subspaces and they would be projected out if a measurement of the syndrome 
is performed. So from the point of view of the basis of the full Hilbert space, the overall transformation under this 
procedure is of the general fault-tolerant type we described. Since the code space is a subspace of an eigenspace of 
the Hamiltonian, it will be effectively transformed by the corresponding transversal operation. 

The standard fault-tolerant procedures provide prescriptions of how to implement any encoded gate by a sequence of 
elementary transversal operations. Therefore, if we make the code space follow an appropriate sequence of transversal 
operations in the described adiabatic manner, when we complete an encoded operation we will have taken the code 
space around a loop whose associated holonomy is equal to the encoded operation. A simple way to see that this 
is indeed a holonomy is to notice that if we track the initial code space as it evolves, it undergoes a loop in the 
Grassmannian since at the end we complete an encoded gate. Furthermore, at all times the Hamiltonian acts trivially 
on the code space so that all states inside it acquire the same dynamical phase. Hence the nontrivial transformation 
resulting inside the code space must be geometric. Thus by following precisely the sequence of transversal operations 
and measurements that are used in a given dynamical fault-tolerant scheme, we obtain a scheme that implements 
logical gates through adiabatic holonomies and at the same time is fault-tolerant. 

In the case of subsystem codes, the code subsystem can be thought of as a collection of subspaces, each of which 
contains the same redundant information. The relative phases between each of these subspaces are gauge degrees of 
freedom. Applying a particular encoded gate is equivalent to applying the same gate in each subspace. In this case, 
our scheme can use Hamiltonians that are either elements of the stabilizer or elements of the gauge group. If the 
Hamiltonian is an element of the stabilizer, all "redundant" subspaces inside a given syndrome subspace belong to a 
single eigenspace of the Hamiltonian, and all of them will undergo the transversal operation that effectively transforms 
that eigenspace. If the Hamiltonian is an element of the gauge group, then some of the subspaces of interest will 
belong to the ground space while others will belong to the excited space. However, since the scheme implements 
the same transversal operation in each eigenspace, all subspaces by construction will undergo the same transversal 
operation. The relative dynamical phases acquired between subspaces in the ground and excited eigenspaces amount 
to a gauge transformation that does not affect the fault tolerance of the operation. Thus by following the same 
sequence of transversal operations and measurements as in a standard dynamical scheme, we obtain a fault-tolerant 
holonomic realization of encoded computation on subsystem codes. 

From the point of view of the full Hilbert space, this method performs transformations of the generalized fault- 
tolerant type we described earlier. This is what ensures the fault tolerance of the method. From a geometrical point 
of view, this corresponds to transporting each of the different syndrome subspaces around a loop such that all of them 
simultaneously undergo the same geometric transformation (in the case of subsystem codes, the same applies for all 
redundant subspaces inside each syndrome space). The statement that all subspaces simultaneously undergo the same 
geometric transformation makes sense with respect to a particular choice of the instantaneous basis in each subspace. 
For any choice of basis in the code space at a given point along the loop, there is a preferred choice of bases inside the 
rest of the syndrome spaces that is determined by the notion of correctable errors. Correctable errors can be thought 
of as transitions from the code space to the other syndrome spaces which can be undone if error correction is applied. 
These errors therefore map the basis from the code space to particular bases in the error spaces such that with respect 
to these bases, a state that has undergone a correctable error has the same form as the form that the non-erroneous 
state has with respect to the basis of the code space. Our scheme applies the same geometric transformation, in this 
sense, in all syndrome subspaces. 

Note that the present approach differs from the original HQC model ™ that it computes in several subspaces 

at the same time. Another difference from the original model is that we do not use a single family of iso-degenerate 
Hamiltonians. This is because for simplicity we use Hamiltonians that are equal to a single element of the stabilizer 
or the gauge group at a given time, and we change the Hamiltonians along different portions of the loop. Thus, if the 
Hamiltonians are elements of the stabilizer, a particular syndrome subspace may belong to the ground space of the 
Hamiltonian during one portion of the loop and to the excited space during another one. Therefore, the holonomies 
in our scheme can be most naturally understood as resulting from parallel transport along loops in the tautological 
bundle over the Grassmannian (this is the fiber bundle whose base is the Grassmannian and whose fibers are the 
subspaces corresponding to the different points in the Grassmannian), rather than loops in a bundle over a space 
parametrizing an iso-degenerate family of Hamiltonians. In the case of subsystem codes, if along the loop we change 
between Hamiltonians which are non-commuting elements of the gauge group, the redundant subspaces that constitute 
a logical subsystem may seem to undergo dynamical transformations in addition to the geometric ones. However, 
these dynamical changes are equivalent to gauge transformations and do not affect the workings of the scheme. We 
could modify our scheme so that it uses a Hamiltonian that separates all of the subspaces of interest by energetic gaps 
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and adiabatically transports each of these subspaces along the same path that it would follow under the scheme we 
described. Then the holonomy resulting in each eigenspace could be understood as being of the original HQC type. 
However, this is unnecessary since the dynamical phases are irrelevant for the workings of our model. 

As single-qubit unitaries together with the C-NOT gate form a universal set of gates, fault-tolerant computation can 
be realized entirely in terms of single-qubit operations and C-NOT operations between qubits from different blocks, 
assuming that the "cat" state can be prepared reliably. Hence our task will be to find adiabatic realizations of these 
operations, as well as of the operations for preparing and verifying the "cat" state. Then these operations can be 
used as building blocks to implement fault-tolerant HQC according to the idea described in this section. 

III. THE SCHEME 

Consider a [[n, l,r, 3]] stabilizer code. This is a code that encodes 1 qubit into n, has r gauge qubits, and can 
correct arbitrary single-qubit errors. To perform a holonomic operation on this code, we need a nontrivial starting 
Hamiltonian that leaves the code space or code subsystem invariant. It is easy to verify that the only Hamiltonians 
that satisfy this property are linear combinations of the elements of the stabilizer and, in the case of subsystem codes, 
elements of the gauge group. 

Note that the stabilizer and the gauge group transform during the course of the computation under the operations 
being applied. At any stage when we complete an encoded operation, they return to their initial forms. Our scheme 
will follow the same transversal operations as those used in a standard dynamical fault-tolerant scheme, but as we 
explained in the previous section, in addition we will have extra dynamical phases that multiply each syndrome 
subspace or are equivalent to gauge transformations. However, it is easy to see that these phases do not affect the 
way the stabilizer or the gauge group transform, so we can omit them from our analysis of the transformation of these 
groups. 

During the implementation of a standard encoded gate, the Pauli group Q„ on a given codeword may change in 
such a way that it acts on other codewords, but it can be verified that this "spreading" can be limited to at most 

4 other codewords including the "cat" state. This is because the encoded C-NOT gate can be implemented fault 
tolerantly on any stabilizer code by a transversal operation on 4 encoded qubits 0, and any encoded Clifford gate 
can be realized using only the encoded C-NOT, provided that we are able to do fault-tolerant measurements (the 
encoded Clifford group is generated by the encoded Hadamard, Phase and C-NOT gates). Encoded gates outside 
of the Clifford group, such as the encoded 7r/8 or Toffoli gates, can be implemented fault tolerantly using encoded 
C-NOT gates conditioned on the qubits in a "cat" state, so they may require transversal operations on a total of 

5 blocks. More precisely, the fault-tolerant implementation of the Toffoli gate requires the preparation of a special 
state of three encoded qubits which involves a sequence of conditional encoded Phase operations and conditional 
encoded C-NOT operations with conditioning on the qubits in a "cat" state Q. But the encoded Phase gate has a 
universal implementation using an encoded C-NOT between the qubit and an ancilla, so the conditional Phase gate 
may require applying a conditional encoded C-NOT. The procedure for implementing an encoded n/8 gate involves 
applying an encoded SX gate conditioned on the qubits in a "cat" state [4^ {S denotes the Phase gate), but the 
encoded S gate generally also involves an encoded C-NOT on the qubit and an ancilla, so it may also require the 
interaction of 4 blocks. For CSS codes, however, the spreading of the Pauli group that acts on a given block can be 
limited to a total of 3 blocks during the implementation of a basic encoded operation since the encoded C-NOT gate 
has a transversal implementation [3]. 

It also should be noted that fault-tolerant encoded Clifford operations can be implemented using only Clifford gates 
on the physical qubits 7\ . These operations transform the stabilizer and the gauge group into subgroups of the Pauli 
group, and their elements remain in the form of tensor products of Pauli matrices. The fault-tolerant implementation 
of encoded gates outside of the Clifford group, however, involves operations that take these groups outside of the 
Pauli group. We will, therefore, consider separately two cases: encoded operations in the Clifford group, and encoded 
operations outside of the Clifford group. 

A. Encoded operations in the Clifford group 

In Ref. 0] it was shown that every encoded operation in the Clifford group can be implemented fault tolerantly 
using Clifford gates on physical qubits. The Clifford group is generated by the Hadamard, Phase and C-NOT gates, 
but in addition to these gates, we will also demonstrate the adiabatic implementation of the X and Z gates which 
are standard for quantum computation. We will restrict our attention to implementing single-qubit unitaries on the 
first qubit in a block, as well as C-NOT operations between the first qubits in two blocks. The operations on the rest 
of the qubits can be obtained analogously. 
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1. Single-qubit unitary operations 



In order to implement a single-qubit operation on the first qubit in a block, we will choose as a starting Hamiltonian 
an element of the stabilizer (with a minus sign) or an element of the gauge group that acts non-trivially on that qubit. 
Since we are considering codes that can correct arbitrary single-qubit errors, one can always find an element G of the 
initial stabilizer or the initial gauge group that has a factor = I , — X , a'^ — Y or = Z acting on the first 
qubit, i.e., 

G = cr*0G, i== 0,1,2,3, (5) 

where G is a tensor product of Pauli matrices and the identity on the remaining n — 1 qubits. It can be verified that 
under Clifford gates the stabilizer and the gauge group transform in such a way that this is always the case except 
that the factor G may spread to qubits in other blocks. From now on, we will use "hat" to denote operators on all 
these qubits and "tilde" to denote operators on all the qubits except the first one. 

Without loss of generality we will assume that the chosen stabilizer or gauge-group element for that qubit has the 
form 

G = Z®G. (6) 

As initial Hamiltonian, we will take the operator 

H{Q)^-G^~Z®G. (7) 

Thus, if G is an element of the stabilizer, the code space will belong to the ground space of H{0). Our goal is to find 
paths in the space of parameters of the Hamiltonian such that when the Hamiltonian is varied adiabatically along 
these paths, each of its eigenspaces undergoes a transformation which is equivalent to that caused by a single-qubit 
operation on the first qubit. 

Proposition 1. If the initial Hamiltonian ([7]) is varied adiabatically so that only the factor acting on the first qubit 
changes, 

H{t) = -H{t) (g) G, (8) 

where 

TT{H{t)} = 0, (9) 



the transformation that each of the eigenspaces of this Hamiltonian undergoes will be equivalent to that resulting 
from a local unitary on the first qubit up to a global phase, i.e., the geometric part Ug{t) = U^^ ® U^_^ of the overall 
unitary, where U\., i = 0,1, are the path-ordered exponents ^ corresponding to the ground and excited spaces, 
respectively, will be equal to a local unitary on the first qubit, Ug{t) = U{t) (g) /. 

Proof. Observe that ([8]) can be written as 

H{t)^H{t)®h~H{t)®A. (10) 

where 

are orthogonal complementary projectors. The evolution driven by H(t) is therefore 

U{t)^Uo{t)(E)Po + Ui{t)(g)Pi, (12) 

where 

t t 

Uo{t)=Texp{-ijH{T)dT), Ui{t) ^Texp{iJ H{T)dT). (13) 
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Let |<?!>o(^)) and \4>i{t)) be the instantaneous ground and excited states of H{t) with eigenvalues Eoit) = —E{t), 
Ei(t) = E{t) {E{t) > 0). Using Eq. ^ for the expressions p^ . we obtain that in the adiabatic hmit, 

C/,(t)=e*'^(*)f/A^.(t)®e-'"Wc/^^(i), J =0,1, (14) 

where Lu{t) — J* dTEij) and 

C/A,W = e/o'^-(^.MI^I^.M)|c/,^.(t))(</>^.(0)|, J -0,1. (15) 
The projectors on the ground and excited eigenspaces of H{Q) are 

Po = |0o(O))(0o(O)|®P^ + |(/'i(O))(0i(O)|®P^ (16) 

and 

Pi - |0i(O))(0i(O)|®Pg + |0„(o))(0o(o)|®P^, (17) 
respectively. Using Eq. and Eq. one can see that the effect of the unitary ([T^ on each of these projectors is 

U{t)Pa^e'^^'\UAo{t)®UA,m®IPo, (18) 

C/(i)Pi =e-*'^(*)(t^A„(i)©C/Ai(t))0 7Pi, (19) 
i.e, up to an overall dynamical phase its effect on each of the eigenspaces is the same as that of the unitary 

Ug{t) =U{t)(g)I, (20) 

where 

U{t) = UAo{t)®UAAt)- (21) 

This completes the proof. 

We next show how by suitably choosing H{t) we can implement all necessary single-qubit gates. We will identify a 
set of points in parameter space, such that by interpolating between these points we can draw various paths resulting 
in the desired transformations. We remark that if a path does not form a loop, the geometric transformation (j2ip 
could be associated with an open-path holonomy taking place inside each eigenspace of the Hamiltonian, provided 
that the final eigenspaces have non-zero overlap with the corresponding initial ones [29| (see Appendix |B| . 

Consider the single-qubit unitary operator 



where 6* is a real parameter. Note that V^^ = {V^^Y . Define the following single-qubit Hamiltonian: 

^e± ^ ye±^y0T. (23) 

Let H{t) in Eq. ([8]) be a Hamiltonian which interpolates between H{0) = Z and H{T) — H^^ (up to a factor) as 
follows: 

H{t) = fit)Z + gm'^^H<l^,it), (24) 

where /(O) > 0, g{T) > 0, and f{T) = g{0) ~ 0. To simplify our notation, we will drop the indices / and g of the 
Hamiltonian, since the exact form of these functions is not important f or our analysi s as long as they are sufficiently 
smooth (see discussion below). This Hamiltonian has eigenvalues ±y^ f{t)^ + g^t)"^ and its energy gap is non-zero 
unless the entire Hamiltonian vanishes. 

Proposition 2. In the adiabatic limit, the Hamiltonian ([8]) with H{t) = H^^{t) gives rise to the geometric transfor- 
mation V^^ ® I. 

The proof of this proposition is presented in Appendix A. 

We will use this result to construct a set of standard gates by sequences of operations of the form V^^ , which can 
be generated by interpolations of the type (j24p run forward or backward. For single-qubit gates in the Clifford group, 
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we will only need three values oi 9: 0, 7r/2 and 7r/4. For completeness, however, we will also demonstrate how to 
implement the tt/S gate, which together with the Hadamard gate is sufficient to generate any single-qubit unitary 
transformation [4^|. For this we will need 6 — n/8. Note that 

H'^^ =±{coseX + smeY), (25) 

so for these values of 9 we have i?o± = ±X, iJ^/2± ^ ±Y, i7'^/4± = ±i^X + -j^Y), H""/^^ = ±(cos f X + sin fF). 
Consider the following adiabatic interpolations: 

-Z®G-^-Y®G-^Z®G. (26) 

According to the above result, the first interpolation yields the transformation The second interpolation can 

be regarded as the inverse oi Z ® G ^ —Y ® G which is equivalent to —Z (g) G Y (g) G since H{t) and —H{t) yield 
the same geometric transformations. Thus the second interpolation results in (l/'^/2^)t — yr/2+ r^Y^^ result is 
therefore ^^72+^^/2+ = jx. We see that up to a global phase, the geometric part of the transformation resulting 
from the above sequence is equal to the single-qubit X gate. 

Similarly, one can verify that the Z gate can be realized via the loop 

-Z®G~^~X®G~^Z®G-^Y®G-^-Z®G. (27) 
The Phase gate can be realized by applying 

-Z®G^-[^X + ^Y)®G^ Z®G, (28) 
V 2 V 2 

followed by the X gate. 

The Hadamard gate can be realized by first applying followed by 

- Z ®G -X ®G. (29) 
Finally, the tt/S gate can be implemented by first applying Y = iXZ, followed by 

Z (gG ^ ~{cos-X + sin-Y)®G ~Z ®G. (30) 
8 8 

Comment. We emphasize that the elementary gates we construct here are not holonomies associated with the code 
space or the error subspaces. The holonomy in a given subspace is a transformation whose components are defined 
with respect to a basis of that subspace, while these elementary gates are the geometric parts of the unitaries resulting 
in the full Hilbert space under the described adiabatic evolutions. As explained in Sec. II C, only after we compose 
a suitable sequence of such elementary gates do we obtain a closed-loop holonomy in the code space (or the code 
subsystem) . The non-Abelian holonomies through which we perform HQC in the code space and the error subspaces 
are the universal set of encoded gates that we would obtain by composing the set of elementary one- and two-qubit 
gates (combined with measurements if necessary) according to the rules of a given dynamical fault-tolerant scheme. 



2. Note on the adiabatic condition 



Before we show how to implement the C-NOT gate between two qubits, let us comment on the conditions under 
which the adiabatic approximation assumed in the above operations is satisfied. Because of the form (|12p of the 
overall unitary, the adiabatic approximation depends on the extent to which each of the unitaries (jl3p approximates 
the expression ([T^ . The latter depends only on the adiabatic properties of the non-degenerate two- level Hamiltonian 
H{t). For such a Hamiltonian, the simple version of the adiabatic condition [4Q| reads 



where 



^ « 1> (31) 



e= max^|(</.i(t)|^|0o(i))l, (32) 
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and 



mm (Eiit) - Eo(t)) = mm 2E(t) (33) 

0<t<T 0<t<T 



is the minimum energy gap of H{t). 

Along the segments of the parameter paths we described, the Hamiltonian is of the form (|24p and its derivative is 

dH^^(t) dfit) doit) , , 

dt dt dt ^ ' 

This derivative is well defined as long as ^^^^ and '^^j^p are well defined. The curves we described, however, may not 
be differentiable at the points connecting two segments. In order for the Hamiltonians that interpolate between 
these points to be differentiable, the functions f{t) and g(t) have to satisfy ^^^p- = and ^^^^ = 0. This means 
that the change of the Hamiltonian slows down to zero at the end of each segment (except for a possible change in 
its strength), and increases again from zero along the next segment. We point out that when the Hamiltonian stops 
changing, we can turn it off completely by decreasing its strength. This can be done arbitrarily fast and it would not 
affect a state which belongs to an eigenspace of the Hamiltonian. Similarly, we can turn on another Hamiltonian for 
the implementation of a different operation. 

The above condition guarantees that the adiabatic approximation is satisfied with precision 0((-^)^). It is known, 
however, that under certain conditions on the Hamiltonian, we can obtain better results (47l . l48j . Let us write the 
Schrodinger equation as 

i|i^(t)) = nmit)) ^ lm\m), m 

where e > is small. Assume that H{t) is smooth and all its derivatives vanish at the end points t = and t ^ T 
(note that this means that H{t) is non-analytic at these points, unless it is constant; an alternative strategy is to 
consider analytic Hamiltonians with a finite number of vanishing derivatives at t = 0,T [48|V Then if we keep H{t) 
fixed and vary e, the adiabatic error would scale super-polynomially with e, i.e., the error will decrease with e faster 
than O(e^) for any N [47]. (Notice that cx e, i.e., the error according to the standard adiabatic approximation is 

In our case, the smoothness condition translates directly to the functions f{t) and g(t). For any smooth f{t) and 
g(t) we can further ensure that the condition at the end points is satisfied by the reparametrization f{t) f{y{t)), 
g{t) g{y{t)) where y{t) is a smooth function of t which satisfies y(0) = 0, y(T) = T, and has vanishing derivatives 
at t = and t = T. Then by slowing down the change of the Hamiltonian by a constant factor e, which amounts to 
increasing the total time T by a factor 1/e, we can decrease the error super-polynomially in e. We will use this result 
to obtain a low-error interpolation in Sec. IV, where we estimate the time needed to implement a geometric gate with 
a certain precision. 



3. The C-NOT gate 

The stabilizer (or gauge group) on multiple blocks of the code is a direct product of the stabilizers (or gauge groups) 
of the individual blocks. Therefore, from Eq. ([5]) it follows that one can always find an element of the initial stabilizer 
or gauge group on multiple blocks that has any desired combination of factors cr*, i = 0, 1,2,3, on the first qubits 
in these blocks. It can be verified that applying transversal Clifford operations on the blocks docs not change this 
property. Therefore, we can find an element of the stabilizer or the gauge group that has the form ([5]), where the 
factor Z acts on the target qubit and G acts trivially on the control qubit. We now explain how to implement the 
C-NOT gate geometrically starting from such a Hamiltonian. 

Notice that a Hamiltonian of the form 

H{t) = \0){0\'' (E)Ha{t)^G+\l){l\'' ^Hi{t)(g)G, (36) 
where the superscript c denotes the control qubit, gives rise to the unitary transformation 

U{t) = \0){0\'®Uoit) + \l){l\'^Ui{t), (37) 

where 

t 

Uasit) = Texp{-i J dTHa,i{T) ® G). (38) 
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If Ho{t) and Hi{t) have the same non-degenerate instantaneous spectra, and Tr{_H'o,i(i)} = 0> then from Eq. ([T]) 

and Proposition 1 it follows that in the adiabatic limit each of the eigenspaces of H{t) will undergo the geometric 
transformation 

Ug{t) ^ |0>(0|"®Vo(i)«)7-l-|l>(l|"®Fi(t)«)7, (39) 

where Vo,i(t) (8) / are the geometric transformations generated by ffo,i(i) ^ G according to Proposition 1. 

Our goal is to find Ho{t) and Hi{t), Hq{0) = Hi{0) = Z, such that at the end of the transformation, the geometric 
unitary ([55)1 will be equal to the C-NOT gate. In other words, we want Vb(2T) — I and Vi{2T) = X (here we have 
chosen the total time of evolution to be 2T for convenience). 

We already saw how to generate geometrically the X gate up to a phase — Eq. (pS)) . We can use the same Hamiltonian 
in place of Hi{t): 

^H''^^-{2T -t), T<t<2T. ' 

Now we want to find a Hamiltonian H(){t) with the same spectrum as Hi(t), which gives rise to a trivial geometric 
transformation, Vo{2T) — I (possibly up to a phase, which can be undone later). Since all Hamiltonians of the type 
H^^{t) have the same instantaneous spectrum (for fixed f{t) and g{t)), we can simply choose 

°^ ' [ff'^/2+(2T-<), T<t<2T, ^ ' 

which corresponds to applying a given transformation from t = to t = T and then undoing it (running it backwards) 
from t = T to t = 2T. This results exactly in Vo{2T) = I. 

Since, as we saw in Sec. III.A.l, the Hamiltonian Hi{t) (g) G gives rise to the geometric transformation iX (g) /, the 
above choice for the Hamiltonians (|41]l and (|40|) in Eq. ([36|l will result in the geometric transformation 

|0)(0|'=(8)/® 7+i|l)(l|'=(gX(g7, (42) 

which is the desired C-NOT gate up to a Phase gate on the control qubit. We can correct the phase by applying the 
inverse of the Phase gate to the control qubit, either before or after the described transformation. 

Notice that from t — Oto t ~ T the Hamiltonians (|40l) and (|4ip are identical, i.e., during this period the Hamiltonian 
(|36l) has the form 

r^H-^/^+{t)(g,G, (43) 

so we are simply applying the single-qubit operation 1/'^/^+ to the target qubit according to the method for single- 
qubit gates described before. It is straightforward to verify that during the second period, from t = T to t = 2T, the 
Hamiltonian (|36p realizes the interpolation 

- r (g) y (g) G -Z" ®Z®G, (44) 

which is understood as in Eq. ([24]) . 

To summarize, the C-NOT gate can be implemented by first applying the inverse of the Phase gate (5^) on the 
control qubit, as well as the transformation on the target qubit, followed by the transformation ([33]). Due to 

the form p6p of H{t), the extent to which the adiabatic approximation is satisfied during this transformation depends 
only on the adiabatic properties of the single-qubit Hamiltonians H^/'^^{t) which we discussed in the previous section. 

Our construction allowed us to prove the resulting geometric transformations without explicitly calculating the 
path-ordered integrals 0. It may be instructive, however, to demonstrate this calculation for at least one of the gates 
we described. In Appendix B, we present an explicit calculation of the geometric transformation for the Z gate for the 
following two cases: f{t) — \ — g(t) = ^ (linear interpolation); f{t) = cos g{t) = sin ^ (unitary interpolation). 

B. Encoded operations outside of the Clifford group 

For universal fault-tolerant computation we also need at least one encoded gate outside of the Clifford group. The 
fault-tolerant implementation of such gates is based on the preparation of a special encoded state [3i H. S EE, 1111 



11 



which involves a measurement of an encoded operator in the Chfford group. For example, the 7r/8 gate requires the 
preparation of the encoded state Io)+°^p(i^/'^)|i) ^ which can be realized by measuring the encoded operator e~™/'^SX 
[4^ . Equivalently, the state can be obtained by applying the encoded operation RS'^ , where R denotes the Hadamard 
gate, on the encoded state |o)-Niin(7r/8) 1 1) .^j^j^,]^ ^^^^ prepared by measuring the encoded Hadamard gate 

0. The Toffoli gate requires the preparation of the three-qubit encoded state looo)+lo^o)+|ioo)+|m) ^nd involves a 
similar procedure [sij . In all these instances, the measurement of the encoded Clifford operator is realized by applying 
transversally the operator conditioned on the qubits in a "cat" state. 

We now show a general method that can be used to implement geometrically any conditional transversal Clifford 
operation with conditioning on the "cat" state. Let O be a Clifford gate acting on the first qubits from some set of 
blocks. As we discussed in the previous section, under this unitary the stabilizer and the gauge group transform in 
such a way that we can always find an element with an arbitrary combination of Pauli matrices on the first qubits. If 
we write this element in the form 

G = Gi®G2,...,n, (45) 

where Gi is a tensor product of Pauli matrices acting on the first qubits from the blocks, and G2,...,n is an operator 
on the rest of the qubits, then applying O conditioned on the first qubit in a "cat" state transforms this stabilizer or 
gauge-group element as follows: 

I'®Gi^G2^...^n = |0)(0|^(8Gi®G2,...,„+|l)(l|"(8Gi(»G2,...,„ |0)(0|=(^Gi(»G2,...,„+|l)(ir®OGiOt^G2,...,„, (46) 

where the superscript c denotes the control qubit from the "cat" state. We can implement this operation by choosing 
the factor Gi to be the same as the one we would use if we wanted to implement the operation O according to the 
previously described procedure. Then we can apply the following Hamiltonian: 

Hc{o)it) = -|0)(0r ®Gi® G2,...,„ - «(t)|l)(ir ® Ho{t) ® G2,...,„, (47) 

where —Ho{t)®G2,...,n is the Hamiltonian that we would use for the implementation of the operation O, and a{t) is a 
real parameter chosen such that at every moment the operator a[t)\l) ® Ho{t)®G2,...,n has the same instantaneous 
spectrum as the operator |0)(0|'^ (g) Gi (g) G2,...,n- This guarantees that the overall Hamiltonian is degenerate and the 
geometric transformation of each of its eigenspaces is given by the operator 

Ug{t) = |0)(0r®/i®/2,...,„ + |l>(ir§5{/o(t)®/2,...,„, (48) 

where Uo{t) is the geometric transformation on the first qubits generated by —Ho{t) ® G2,...,n- Since we presented 
the constructions of our basic Clifford operations up to an overall phase, the operation Uo{t) may differ from the 
desired operation by a phase. This phase can be corrected by applying a suitable gate on the control qubit from the 
"cat" state (we explain how this can be done in the next section). We remark that a Hamiltonian of the type ((47|) 
requires fine tuning of the parameter a{t) and generally can be complicated. Our goal in this section is to prove that 
universal fault-tolerant holonomic computation is possible in principle. In Sec. V we show that depending on the 
code one can find more natural implementations of these operations. 

If we want to apply a second conditional Clifford operation Q on the first qubits in the blocks, we can do this as 
follows. Imagine that if we had to apply the operation Q following the operation O, we would use the Hamiltonian 
HQ{t) — —HQ{t) ® G2. „, where Hq{0) — OG'iO^ (Xi G2^ „ is a suitable element of the stabilizer or the gauge group 
after the application of O. Before the application of O, that element would have had the form G'l ^G'2 „. Under the 
application of a conditional O, the element Gi(g)G2^ „ transforms to |0)(0|'=(g)G'i(g)G2_ _„-l-|l)(l|'=(g)OGiO'f(g)G2_ „ 
which can be used (with a minus sign) as a starting Hamiltonian for a subsequent operation. In particular, we can 
implement the conditional Q following the conditional O using the Hamiltonian 

HciQ){t) = -|0)(0r ® G'l ® G'2_„ - /3(i)|l)(ir ® HQ{t) ® G'2_^, (49) 

where the factor P{t) guarantees that there is no splitting of the energy levels. Subsequent operations can be applied 
analogously. Using this general method, we can implement a unitary whose geometric part is equal to any transversal 
Clifford operation conditioned on the "cat" state. 

C. Preparing and using the "cat" state 

In addition to transversal operations, a complete fault-tolerant scheme requires the ability to prepare, verify and 
use a special ancillary state such as the "cat" state (|00...0) -I- |ll...l))/\/2 proposed by Shor 2]. This can also be 
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done using our geometric approach. Since the "cat" state is known and its construction is non- fault-tolerant, we 
can prepare it by simply treating each initially prepared qubit as a simple code (with G in Eq. ^ being trivial), 
and updating the stabilizer of the code via the applied geometric transformation as the operation progresses. The 
stabilizer of the prepared "cat" state is generated by ZiZj, i < j. Transversal unitary operations between the "cat" 
state and other codewords are applied as described in the previous section. 

We also have to be able to measure the parity of the state, which requires the ability to apply successive C-NOT 
operations from two different qubits in the "cat" state to the same ancillary qubit initially prepared in the state |0). 
We can regard a qubit in state |0) as a simple code with stabilizer (Z), and we can apply the first C-NOT as described 
before. Even though after this operation the state of the target qubit is unknown, the second C-NOT gate can be 
applied via the same interaction, since the transformation undergone by each eigenspace would still be equivalent to 
the desired C-NOT and at the end when we measure the qubit we project onto one of the eigenspaces. 

D. Fault tolerance of the scheme 

We showed how we can generate any transversal operation on the code space geometrically, assuming that the state 
is non-erroneous. But what if an error occurs on one of the qubits? 

At any moment, we can distinguish two types of errors — those that result in transitions between the ground and 
the excited spaces of the current Hamiltonian, and those that result in transformations inside the eigenspaces. Due 
to the discretization of errors in QEC, it suffices to prove correctability for each type separately. The key property of 
our construction is that the transformation undergone by each of the eigenspaces is equivalent to the same transversal 
operation. Because of this, if we are applying a unitary on the first qubit, an error on that qubit will remain localized 
regardless of whether it causes an excitation or not. If the error occurs on one of the other qubits, at the end of the 
transformation the result would be the desired single-qubit unitary gate plus the error on the other qubit, which is 
correctable. 

It is remarkable that even though the Hamiltonian couples qubits within the same block, single-qubit errors do not 
propagate. This is because the coupling between the qubits amounts to a change in the relative phase between the 
ground and excited spaces, but the latter is irrelevant since either it is equivalent to a gauge transformation, or when 
we apply a correcting operation we project onto one of the eigenspaces. In the case of the C-NOT gate, an error can 
propagate between the control and the target qubits, but it never results in two errors within the same codeword. 

IV. EFFECTS ON THE ACCURACY THRESHOLD FOR ENVIRONMENTAL NOISE 

Since the method we presented conforms completely to a given fault-tolerant scheme, it would not affect the error 
threshold per operation for that scheme. Some of its features, however, would affect the threshold for environment 
noise. 

First, observe that when applying the Hamiltonian ([H]), we cannot at the same time apply operations on the other 
qubits on which the factor G acts non-trivially. Thus, some operations at the lowest level of concatenation that would 
otherwise be implemented simultaneously might have to be implemented serially. The effect of this is equivalent to 
slowing down the circuit by a constant factor. (Note that we could also vary the factor G simultaneously with H{t), 
but in order to obtain the same precision as that we would achieve by a serial implementation, we would have to slow 
down the change of the Hamiltonian by the same factor.) The slowdown factor resulting from this loss of parallelism 
is usually small since this problem occurs only at the lowest level of concatenation. For example, when implementing 
encoded single-qubit gates with the Bacon-Shor code, we can apply operations on up to 6 out of the 9 qubits in a 
block simultaneously. As we show in Sec. V, we can address any two qubits in a row or column using our method 
by taking G in Eq. ([5]) to be a single-qubit operator Z or X on the third qubit in the same row or column. The 
Hamiltonians used to apply operations on the two qubits commute with each other at all times and do not interfere. 
A similar phenomenon holds for the implementation of the encoded C-NOT gate, or the operations involving the 
"cat" state. Thus, for the Bacon-Shor code we have a slowdown due to parallelism by a factor of 1.5. 

A more significant slowdown results from the fact that the evolution is adiabatic. In order to obtain a rough 
estimate of the slowdown due specifically to the adiabatic requirement, we will compare the time Th needed for the 
adiabatic implementation of a given gate with precision 1 — d to the time Td needed for a dynamical realization of the 
same gate with the same strength of the Hamiltonian. We will consider a realization of the X gate via the unitary 
interpolation [ioj 

H{t) = -VxiT{t))ZV^{T{t))<i^G, Vx{Tit))^exp(tT{t)^x] , (50) 
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where r(0) = 0, T{Th) = Th- Thus the energy gap of the Hamiltonian is constant. The optimal dynamical implemen- 
tation of the same gate is via the Hamiltonian —X for time Tj, — ^. 

As we argued in Sec. Ill, the accuracy with which the adiabatic approximation holds for the Hamiltonian ([50|) is 
the same as that for the Hamiltonian 

H{t)=Vx{Tit))ZVUr{t)). (51) 
We now present estimates for two different choices of the function T{t). The first one is 

T{t)^t. (52) 



In this case the Schrodinger equation can be easily solved in the instantaneous eigenbasis of the Hamiltonian ([51 
For the probability that the initial ground state remains a ground state at the end of the evolution, we obtain 



1 

P = T— 



where 



_^cos2(Jv/TT7^) = l-5, (53) 

e = Y^- (54) 

Expanding in powers of e and averaging the square of the cosine whose period is much smaller than Th, we obtain 
the condition 

< 25. (55) 

Assuming, for example, that 5 « 10^** (approximately the threshold for the 9-qubit Bacon-Shor code [i3|)) we obtain 
that the time of evolution for the adiabatic case must be about 70 times longer than that in the dynamical case. 

It is known, however, that if H{t) is smooth and all its derivatives vanish at i = and t ^ Th, the adiabatic error 
decreases super-polynomially with Th [47,]. To achieve this, we will choose 

T(t) = i /'dt'e-i/^'^^"*'/^''), dt'e-^^^'"'^^''/^^\ (56) 



For this interpolation, by a numerical solution we obtain that when Th/Td ~ 17 the error is already of the order 
of 10^®, which is well below the threshold values obtained for the Bacon-Shor codes This is a remarkable 

improvement in comparison to the previous interpolation which shows that the smoothness of the Hamiltonian plays 
an important role in the performance of the scheme. 

An additional slowdown in comparison to a perfect dynamical scheme may result from the fact that the constructions 
for some of the standard gates we presented involve long sequences of loops. With more efficient parameter paths, 
however, it should be possible to reduce this slowdown to minimum. An approach for finding optimal loops presented 
in Ref. [50| niay be useful in this respect. 

In comparison to a dynamical implementation, the allowed rate of environmental noise for the holonomic case would 
decrease by a factor similar to the slowdown factor. In practice, however, dynamical gates are not perfect and the 
holonomic approach may be advantageous if it gives rise to a better operational precision. 

We finally point out that an error in the factor H{t) in the Hamiltonian ([5]) would result in an error on the first 
qubit according to Eq. (|2T|) . Such an error clearly has to be below the accuracy threshold. More dangerous errors, 
however, are also possible. For example, if the degeneracy of the Hamiltonian is broken, this can result in an unwanted 
dynamical transformation affecting all qubits on which the Hamiltonian acts non-trivially. Such multi-qubit errors 
have to be of higher order in the threshold, which imposes more severe restrictions on the Hamiltonian. 



V. FAULT-TOLERANT HOLONOMIC COMPUTATION WITH LOW- WEIGHT HAMILTONIANS 

The weight of the Hamiltonians needed for the scheme we described depends on the weight of the stabilizer or 
gauge-group elements. Remarkably, certain codes possess stabilizer or gauge-group elements of low weight covering 
all qubits in the code, which allows us to perform holonomic computation usin g lo w-weight Hamiltonians. Here we 
will consider as an example a subsystem generalization of the 9-qubit Shor cod e fiol — the Bacon-Shor code [2^[23] — 
which has particularly favorable properties for fault-tolerant computation [40l |52{ . In the 9-qubit Bacon-Shor code, 
the gauge group is generated by the weight-two operators Zk.jZk,j+i and Xj,kXj^i^k, where the subscripts label the 
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qubits by row and column when they are arranged in a 3 x 3 square lattice. Since the Bacon-Shor code is a CSS code, 
the C-NOT gate has a direct transversal implementation. We now show that the C-NOT gate can be realized using 
at most weight-three Hamiltonians. 

If we want to apply a C-NOT gate between two qubits each of which is, say, in the first row and column of its 
block, we can use as a starting Hamiltonian —Z\ ^ ^ ZI2, where the superscript t signifies that these are operators 
in the target block. We can then apply the C-NOT gate as described in Sec. III. After the operation, however, this 
gauge-group element will transform io — Z\ ^ ® ^ ® Z\ 2- If we now want to implement a C-NOT gate between 
the qubits with index {1, 2} using as a starting Hamiltonian the operator —Z\ Z^-^® Z\2 according to the same 
procedure, we will have to use a four-qubit Hamiltonian. Of course, at this point we can use the starting Hamiltonian 
—Z\ 2 ^ Zl 2, but if we had also applied a C-NOT between the qubits labeled {1,3}, this operator would not be 
available — it would have transformed to — Z\ 2 ® Z\ ® Zl 

What we can do instead, is to use as a starting Hamiltonian the operator — i ^ Zl 2 'S) Zf 2 which is obtained 
from the gauge-group element Zl ^ (g) Zf i (g) 2 ® ■Z'1.2 after the apphcation of the C-NOT between the qubits with 

index {1, 1}. Since the C-NOT gate is its own inverse, we can regard the factor Z* as G in Eq. and use this 
starting Hamiltonian to apply our procedure backwards. Thus we can implement any transversal C-NOT gate using 
at most weight-three Hamiltonians. 

Since the encoded X, Y and Z operations have a bitwise implementation, we can always apply them according to 
our procedure using Hamiltonians of weight 2. For the Bacon-Shor code, the encoded Hadamard gate can be applied 
via bitwise Hadamard transformations followed by a rotation of the grid by a 90 degree angle [i^l . The encoded Phase 
gate can be implemented by using the encoded C-NOT and an ancilla. 

The preparation and measurement of the "cat" state can also be done using Hamiltonians of weight 2. To prepare 
the "cat" state, we first prepare all qubits in the state (|0) -I- |l))/-\/2, which can be done by measuring each of 
them in the {|0),|1)} basis (this ability is assumed for any type of computation) and applying the transformation 
— Z —X or Z ^ —X depending on the outcome. To complete the preparation of the "cat" state, apply a two-qubit 
transformation between the first qubit and each of the other qubits {j > 1) via the transformation 

-h^Xj-^ -Zi(E)Zj. (57) 

Single-qubit transformations on qubits from the "cat" state can be applied according to the method described in the 
previous section using at most weight-two Hamiltonians. 

To measure the parity of the state, we need to apply successively C-NOT operations from two different qubits in 
the "cat" state to the same ancillary qubit initially prepared in the state |0). As described in Sec. HI, this can also 
be done according to our method and requires Hamiltonians of weight 2. 

For universal computation with the Bacon-Shor code, we also need to be able to apply one encoded transformation 
outside of the Clifford group. As we mentioned earlier, in order to implement the Toffoli gate or the tt/8 gate, it is 
sufficient to be able to implement a C-NOT gate conditioned on a "cat" state. For the Bacon-Shor code, the C-NOT 
gate has a transversal implementation, so the conditioned C-NOT gate can be realized by a series of transversal Toffoli 
operations between the "cat" state and the two encoded states. We now show that this gate can be implemented 
using at most three-qubit Hamiltonians. 

Ref. [33I provides a circuit for implementing the Toffoli gate as a sequence of one- and two-qubit gates. We will use 
the same circuit, except that we flip the control and target qubits in every C-NOT gate using the identity 

(i?l «)i?2)Ci,2(i?l ® i?2) = C2,l, (58) 

where Ri denotes a Hadamard gate on the qubit labeled by i and Cij denotes a C-NOT gate between qubits i and 
j with i being the control and j being the target. Let Toffoli^ ^ ^ denote the Toffoli gate on qubits i, j and k with i 
and j being the two control qubits and k being the target qubit, and let Si and Ti denote the Phase and 7r/8 gates 
on qubit i, respectively. Then the Toffoli gate on three qubits (the first one of which we will assume to belong to the 
"cat" state), can be written as: 

Toffolii,2,3 = R2Cz^2RzTlR^RiCz,iR3T^R?,C3^2R?,TlR-iCz.iR^T-iR^R2TlR2C2,iR2TlR2C2,iR2S2RiTi. (59) 

To show that each of the above gates can be implemented according to our geometric approach using Hamiltonians 
of weight at most 3, we will need an implementation of the C-NOT gate which is suitable for the case when we have 
a stabilizer or gauge-group element of the form 



G = X®G, 



(60) 
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where the factor X acts on the target qubit and G acts trivially on the control qubit. By a similar argument to the 
one in Sec. Ill, one can verify that in this case the C-NOT gate can be implemented as follows: apply the operation 
on the control qubit (we describe how to do this for our particular case below) together with the transformation 

-X®G^-Z®G^X®G (61) 

on the target qubit, followed by the transformation 

r «) X (g) G ^ -(|0)(0|" ®Z+\l){l\''®Y)®G^-F®X®G. (62) 

Since the second and the third qubits belong to blocks encoded with the Bacon-Shor code, there are weight-two 
elements of the initial gauge group of the form Z ® Z covering all qubits. The stabilizer generators on the "cat" state 
are also of this type. Following the transformation of these operators according to the sequence of operations (j59p . 
one can see that before every C-NOT gate in this sequence, there is an element of the form ([TO]) with G = Z that 
can be used to implement the C-NOT gate as described, provided that we can implement the gate on the control 
qubit. We also point out that all single-qubit operations on qubit 1 in this sequence can be implemented according 
to the procedure described in Sec. Ill, since at every step we have a weight-two stabilizer element on that qubit with 
a suitable form. Therefore, all we need to show is how to implement the necessary single-qubit operations on qubits 
2 and 3. Due to the complicated transformation of the gauge-group elements during the sequence of operations (j59p . 
we introduce a method of applying a single-qubit operation with a starting Hamiltonian that acts trivially on the 
qubit. For implementing single-qubit operations on qubits 2 and 3 we use as a starting Hamiltonian the operator 

H{Q) = -Ii®Xi®Z, i = 2,3, (63) 

where the first factor (7^) acts on the qubit on which we want to apply the operation (2 or 3), and Xi ® Z is the 
transformed (after the Hadamard gate i?i) stabilizer element of the "cat" state that acts non-trivially on qubit 1 (the 
factor Z acts on some other qubit in the "cat" state). 

To implement a single-qubit gate on qubit 3 for example, we first apply the interpolation 

-h®Xi®Z^ -Z-i ®Zi®Z. (64) 

This results in a two-qubit geometric transformation on qubits 1 and 3. We do not have to calculate this 
transformation exactly since we will undo it later, but the fact that each eigenspace undergoes the same two-qubit 
geometric transformation can be verified similarly to the C-NOT gate we described in Sec. III. 

At this point, the Hamiltonian is of the form ([7|) with respect to qubit 3, and we can apply any single-qubit unitary 
gate V3 according to the method described in Sec. III. This transforms the Hamiltonian to —V^Z^V^ (S) Zi (S) Z. We 
can now "undo" the transformation C/1.3 by the interpolation 

- V3Z3V^ ® Zi®Z ^ -h® Xi®Z. (65) 

The latter transformation is the inverse of Eq. (|64p up to the single-qubit unitary transformation V3, i.e., it results in 
the transformation V^Ul^V^. Thus the net result is 

VzUl^V^ViUi,^ = ^3, (66) 

which is the desired single-qubit unitary transformation on qubit 3. We point out that during this transformation, 
a single-qubit error can propagate between qubits 1 and 3, but this is not a problem since we are implementing a 
transversal Toffoli operation and such an error would not result in more than one error per block of the code. 

We showed that for the Bacon-Shor code our scheme can be implemented with at most 3-local Hamiltonians. 
This is optimal for the construction we presented, since there are no non-trivial codes with stabilizer or gauge-group 
elements of weight smaller than 2 covering all qubits. One could argue that since the only Hamiltonians that leave 
the code space invariant are superpositions of elements of the stabilizer or the gauge group, one cannot do better than 
this. However, it may be possible to approximate the necessary Hamiltonians with sufficient precision usin g 2- local 
interactions. A possible direction to consider in this respect are the gadget techniques introduced in Ref. [53| and 
developed further in Refs. [13, H^- This is left as a problem for future investigation. 
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VI. CONCLUSION 



We described a scheme for fault-tolerant holonomic computation on stabilizer codes, which demonstrates that HQC 
is a scalable method of computation. The scheme opens the possibility for combining the software protection of error 
correction with the inherent robustness of HQC against control imperfections. Our construction uses Hamiltonians 
that are elements of the stabilizer or the gauge group of the code and works by adiabatically varying these Hamiltonians 
in a manner which generates unitaries whose geometric parts are equal to transversal operations from the point of view 
of the basis of the full Hilbert space. By composing these transversal operations as in a given standard dynamical fault- 
tolerant scheme, we thus transport the code space (or the redundant subspaces that constitute a code subsystem) 
adiabatically around loops that give rise to holonomies equal to encoded operations inside the code space. The 
Hamiltonians needed for implementing two-qubit gates are at least 3-local. We showed that computation with at 
most 3-local Hamiltonians is possible with the Bacon-Shor code. 

It is interesting to point out that the adiabatic regime in which our scheme operates is consistent with the model of 
Markovian decoherence. In Ref. [56| it was argued that the standard dynamical paradigm of fault tolerance is based on 
assumptions that are in conflict with the rigorous derivation of the Markovian limit. Although the threshold theorem 
has been extended to non-Markovian models (stI . [ssI . [s^ , the Markovian assumption is an accurate approximation for 
a wide range of physical scenarios [60| . It also allows for a much simpler description of the evolution in comparison 
to non-Markovian models (see, e.g., Ref. [6l[). In Ref. [5^ it was shown that the weak-coupling-limit derivation of 
the Markovian approximation is consistent with computational methods that employ slow transformations, such as 
adiabatic quantum computation [i^l or HQC. A theory of fault tolerance for the adiabatic model of computation at 
present is not known, although some steps in this direction have been undertaken [H, H^l ■ Our hybrid HQC-QEC 
scheme provides a solution for the case of HQC. We point out, however, that it is an open problem whether the 
Markovian approximation makes sense for a fixed value of the adiabatic slowness parameter when the circuit increases 
in size. Giving a definitive answer to this question requires a rigorous analysis of the accumulation of non-Markovian 
errors due to deviation from perfect adiabaticity. 

Applying the present strategy to actual physical systems might require modifying our abstract construction in 
accordance with the available interactions, possibly using linear combinations of stabilizer or gauge-group elements 
rather than single elements as the basic Hamiltonians. Given that simple QEC codes and two-qubit geometric 
transformations have been realized using NMR (65l.[6^ and ion-trap [67l[68| techniques, these systems seem particularly 
suitable for hybrid HQC-QEC implementations. 

We hope that the techniques presented in this paper might prove useful in other areas as well. It is possible 
that some combination of transversal adiabatic transformations and active correction could provide a solution to the 
problem of fault tolerance in the adiabatic model of computation as well. 
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Appendix A: PROOF OF PROPOSITION 2 

To prove Proposition 2, observe that the Hermitian unitary matrix 

ie- 
-ie'^ 



W = ( '"n ) (Al) 



has the properties 

[ly^y^i] = 0, (A2) 

{W',Z}^0, (A3) 
where [•, •] and {•, •} denote a commutator and an anticommutator, respectively. This means that 

W'^H^^{t)W^ = W\f{t)Z + g{t)V''^ZV^^)W'^ = ~fit)Z - git)V^^ZV'^^ = -^"^(0, (A4) 
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i.e., the unitary flips the ground and excited spaces of H^^{t) for any t. 

The unitaries C/p^, given by Eq. (fT^ for H{t) — H^^{t), are therefore related by 

[/»±(t) = w'^U^,^{t)W^. (A5) 

Using the fact that 

M/^|0) = -ie'^ll), 

W^\l) ^ie"'^\0), (A6) 
from Eq. (I14p and Eq. one can see that Eq. (jASp implies 

[/l±W = W^'f/ltWW^^ (A7) 

Let us define the eigenstates of H''^{t) at time T as |0o'^(r)) = V'^^\0) and |</>i=^(T)) = V^^±|l). Expression ^ can 
then be written as 

[;«±(r) = e-oV«±|0)(0|, 

C/I±(r) = e-'V'^±|1>(1|, (A8) 

where and a^^ are geometric phases. Without explicitly calculating the geometric phases, from Eq. (|A8p . 
Eq. jSI, Eq. jM]) and Eq. (|XT|), we obtain 

e'"o*-l/«±|0)(0| = e*"i*iy«y^±|l)(l|W^^ = e*"i*y^±Ty^|l)(l|W^^ = e*"i*y^±|0)(0|, (A9) 

i.e., 

e'"« =e*"i . (AlO) 

Therefore, up to a global phase, Eq. ([2T|) yields 

L/e±(T)~y«±. (All) 

Appendix B: CALCULATING THE HOLONOMY FOR THE Z GATE 

1. Linear interpolation 

We first demonstrate how to calculate the ground-space holonomy for the Z gate for the case of linear interpolations 
along the segments of the path, i.e., when /(t) and g{t) in Eq. p4)) are 

/W = i-|, .gW-f (Bi) 

In order to calculate the holonomy ([¥]) corresponding to our construction of the Z gate, we need to define a single- 
valued orthonormal basis of the ground space of the Hamiltonian along the loop described by Eq. (j27p . Since the 
Hamiltonian has the form at all times, it is convenient to choose the basis of the form 

\jk;X)^\x,{X)Mk), (B2) 
J = 0,1; fc = l,...,2"-2 , 

where \xo{X{t))) and IxiiM^))) are ground and excited states of H{t), and \ipok) and \ipik) are fixed orthonormal bases 
of the subspaces that support the projectors Pq and Pi defined in Eq. pT|) . respectively. The eigenstates |xo(A(t))) 
and |xi(A(t))) are defined up to an overall phase, but we have to choose the phase such that the states are single- valued 
along the loop. 

Observe that because of this choice of basis, the matrix elements (jS]) become 

(A.W^o' = (jfc;A|^|/fc';A) = {x,{X)\^\xrW) 

~ ~ d 
x(V'ifclV'i'fc') = (Xj(A)|^|Xj'(A))%'4fc', (B3) 
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i.e., the matrix is diagonal. (Since we are looking only at the ground space, we are not writing the index of the 
energy level.) We can therefore drop the path-ordering operator. The resulting unitary matrix C^Jj, j/j,/ acting on the 
subspace spanned by {\jk; A(0))} is also diagonal and its diagonal elements are 



Ujk.k = exp (/(x,(A)|^|x,(A)>dA^) 



(B4) 



These are precisely the Berry phases for the loops described by the states |xi(A)). Since the loop in parameter space 
consists of four line segments, we can write the last expression as 

Ujk,, = exp fe(A)|^|x,(A))dA-^ , (B5) 



where 7^, z = 1, 2, 3, 4, are the segments indexed by their order corresponding to Eq. (j27p . If we parametrize each line 
segment by the dimensionless time < s < 1, we obtain 

U],,,, = exp j\^^^{s)\^^\x){s))d^ , (B6) 
where the superscript i in \x]{s)) indicates the segment. In the {|0), |1)} basis, we will write these states as 

IxK^)) = fSSV J =0,1 , z = 1,2,3,4 , (B7) 



h]{s) 

where |a}(s)p + |6j(s)|2 = 1. 

Along the segment 71, the Hamiltonian has the form H{s) = Hi{s) ® Pq — Hi{s) (g) Pi, where 

Hi{t) = {1- s)Z + sX, (B8) 

i.e, the states |xo(s)) and |x}(s)) are the ground and excited states of i?i(s). For these states we obtain 

1. X (l-. + x/l-2. + 2,s2)e-o(^) 

ao(s) = = , (B9) 

^2-4:3 + 4s2 + (2 - 2s)Vl -2s + 2s2 
blis) = =, (BIO) 





-4s + 4s2 + (2-2s)\/l- 


2s + 2s2 










-4s + 4s2 + (2-2.s)Vl- 


2s + 2s2 




(1-s- Vl-2s + 2s2)e*'^J('*) 




-4s + 4s2- (2-2s)Vl- 


2s + 2s2 










-4s + 4s2- (2-2s)Vl- 


2s + 2s2 



, ' ' (Bll) 

^2 - 4s + 4s2 - (2 - 2s)Vl-2s + 2s2 

blis) = (B12) 
y^2 - 4s + 4s2 - (2 - 2s)Vl-2s + 2s2 

where '^j(s) are arbitrary phases which have to be chosen so that when we complete the loop, the phases of the 
corresponding states will return to their initial values modulo 27r. We will define the loops as interpolating between 
the following intermediate states defined with their overall phases: 

IV'o(A)) : |0) - 1/^) |1) \fZ^') - |0), (B13) 

IVi(A)) : |1) ^ |/°) ^ |0> ^ I/;/') ^ |1), (B14) 

where 

Ui) = (B15) 

In other words, we impose the conditions |xj,i(0)> = !«' !>' IXo,i(l)) = l/±) = IXo.iW), |Xo,i(l)> = |1,0) = IxlM), 

1x^,1(1)) = 1/^^') = 1x^1(0)), lx^,i(i)) = |0'i)- 
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From Eq. (|B9| and Eq. (|BT0)) we see that aj(0) = e^'^o^"', bl{0) = and aj(l) = -i=e*"o(i\ fej(l) = -^e''^"''^^ , so 
we can choose 

wi(s) =0, VsG [0,1]. (B16) 

Similarly, from Eq. ((BTTI) and Eq. ((BT2| it can be seen that a}(0) = 0, 6}(0) = e*"J(") and aj(l) = --i=e''^i'(i\ 

&i(l) — -i=e'"i(^^ This means that ujI{s) has to satisfy e^'^i^"^ = 1, e*"i^^^ = -1. We can choose any differentiable 
i^l{s) that satisfies 

ujI{0)^0, ujI{1) = 7:. (B17) 
In order to calculate (x](s)|^|x](s))rfs, we also need 



\x]{s)) ^ { ^iX'i) . (B18) 



Differentiating Eqs. ((B9|1 - (|BT2)1 yields 



—ar.(s) = = ^ , (B19 

' 2V2 - 4s + 4s2[l -2s + 2s^ + (1 - s)Vl - 2s + 2s^]i 

d_^ = 2-4s + 3s2 + (2-2s)Vl-2s + 2s2 

°^ ^ 2V2 - 4s + 4s2[l - 2s + 2s2 + (1 - s)Vl - 2s + 2s2] I ' 

d 1.. s(l-s-Vl-2s + 2s2)e'"'(^) i d i 

ds ' 2 V2 - 4s + 4s2 [1 - 2s + 2s2 - (1 - s) Vl - 2s + 2s2] i ' ds 

rf^i.^ (2-4s + 3s2-(2-2s)VT^27T2i^)e-i(-) ^ ^^if ( ^32) 

rfs ° 2V2 - 4s + 4s2[l - 2s + 2s2 - (1 - s)Vl - 2s + 2s2]i ^^■')' ds'^^^'''- 



By a straightforward substitution, we obtain 



{xUs)\^Jxl{s)) = aj*(s)^aj(s)+6i*(s)^6j(s) = 0, (B23) 
{xl{s)\-^Jxli^)) = al*is)^^alis)+bl*is)-^^blis)^^-^y,is). (B24) 



Thus the integrals are 







\xlis)\^yois))ds = 0, (B25) 

\xl{s)\^Jxlis))ds = ^u;l{s)\l = ^n. (B26) 
In the same manner, we calculate the contributions of the other three line segments. The results are: 

\xlis)\4:\xlis))ds = 0, (B27) 







ds 

\xlis)\^Jxlim^ = 0, (B28) 



,3/„M \,,3/ 



(Xo(s)|^IXo(s))ds = (B29) 
\xlis)\-^Jxlis))ds = 0, (B30) 
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{Xtis)\^jxti^))ds = 0, (B31) 
' {xtis)\^Jxtis))ds = (B32) 
Putting everything together, for the diagonal elements of the holonomy we obtain 

r/7 _ if 

= (B33) 
The holonomy transforms any state in the ground space of the initial Hamiltonian as 

[/^5]a,fe|j)|V;,,) = e^t 3 = 0,1. (B34) 

jk jk 

From the point of view of the full Hilbert space, this is effectively a Z gate on the first qubit up to an overall phase. 

Note that other single-qubit transformations such as the Hadamard or the X gates, under which the eigenspaces 
of the Hamiltonian do not follow complete loops, can be obtained in a similar fashion by calculating the open-path 
expression In principle, the result of that calculation depends on the choice of basis {\a; A)} which is defined up 
to a unitary gauge transformation. However, if the final eigenspace has a non-zero overlap with the initial one, this 
ambiguity can be removed by defining the frame in the final eigenspace to be the one which is "most parallel" to the 
initial frame poj . In the case of the single-qubit Hadamard gate, there is a non-zero overlap between the initial and 
final eigenspaces, and the set of initial basis states |0A;;0) = |O)|'0ofe), \^k;0) = k — 1,...,2"^^, can be seen 

to be most parallel to the set of final basis states |Ofc; 1) = |-l-)|'0ofe), 1) = |— k = 1, ...,2"~^, respectively 
(for a precise definition of "most parallel" see Ref. [2^). Thus, the resulting open-path holonomy corresponds to 
flipping the phase of half of the basis vectors. For the single-qubit X gate, however, the final ground (excited) space is 
orthogonal to the initial ground (excited) space. Thus a gauge invariant expression for the geometric transformation 
cannot be defined in this way. 

We emphasize that the expression ([2]) for the open-path geometric transformation taking place inside each eigenspace 
of the Hamiltonian, whether gauge-invariant or not, is not the same as the geometric transformation being realized 
inside the code space which is what eventually gives rise to the logical transformation. For example, the single-qubit 
Z gate can be understood as a closed-loop holonomy in the eigenspaces of the Hamiltonian —Z ® G, but under this 
transformation the code space does not follow a loop; in fact, it becomes orthogonal to the initial code space so that 
a gauge-invariant expression for the geometric transformation taking place inside it cannot be defined. However, if a 
sequence of single-qubit Z gates on all qubits yields, say, an encoded Z gate (as in the case of CSS codes), then after 
we complete such a sequence we will obtain a non-trivial closed-loop holonomy in the code space which is equal to 
the encoded Z. 

2. Unitary interpolation 

The calculation is simpler if we choose a unitary interpolation, 

7rt Trt 
/W=cos— , .g(t)=sin— . (B35) 

Such an interpolation corresponds to a rotation of the Bloch sphere around a particular axis for each of the segments 
of the loop. The first two segments of the loop (p7|) are realized via the Hamiltonian 

HiAt) = -VUt)ZVY{t)®G, VV(i) =exp(zt^r) , (B36) 

applied for time T, and the third and fourth segments are realized via the Hamiltonian 

H3A{t) = -Vxit)ZVl{t)(^G, Vx{t)=e^p{it^xy (B37) 

again applied for time T. Let us define the eigenstates of the Hamiltonian along the first two segments as 

\xl''it))^Vx{tm, \x{-'{t))=Vxm), 0<t<T, (B38) 
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and along the third and fourth segments as 

\xl-\t)) = -tVUt)Y\0), \xlHt))=-tVUt)Y\l), 0<t<T. (B39) 

Notice that 

\xl''iT)) = -zY\0) = |x^''(0)), \xl-'iT)) = ~^Y\1) = |x?-'(0)), (B40) 

but 

IXo 'm = |0) ^ \xo Hn - -^|0), |x}''(0)> = |1> ^ \xl'\T)) (B41) 
i.e., this basis is not single- valued. To make it single valued, we can modify it along the third and fourth segments as 
\Xo\t)) ^ \xrit)) = e-o^lxg-'W), \xl-\t)) ^ Ixl'Ht)) - e-^W|xr(i)>, (B42) 

where 

wo(0) = 0, c^o(r) = |, 
o;i(0) = 0, c.i(r) = -| 

The expression (|B6[) then becomes 

f/^.,fc=exp^''(xj''(t)|||xj''W)dt + ^^(xf W>dt + (-lF^ , J =0,1. (B45) 

But 

(X/'W|||X/'W) = -^^Ol^b) = 0, (B46) 

and 

{xf'{t)\j^\xf'it)) = = 0. (B47) 

Therefore, we obtain Eq. (jB33|) . 



(B43) 
(B44) 
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